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The quantitative contributions of a mixed phase-space to 
the mean characterizing the distribution of diagonal transi- 
tion matrix elements and to the variance characterizing the 
distributions of non-diagonal transition matrix elements are 
studied. It is shown that the mean can be expressed as the 
sum of suitably weighted classical averages along an ergodic 
trajectory and along the stable periodic orbits. Similarly, it 
is shown that the values of the variance are well reproduced 
by the sum of the suitably weighted Fourier transforms of 
classical autocorrelation functions along an ergodic trajectory 
and along the stable periodic orbits. The illustrative numeri- 
cal computations are done in the framework of the Hydrogen 
atom in a strong magnetic field, for three different values of 
the scaled energy. 

PACS numbers: 05.45. +b, 03.65.Sq, 32.60. +i, 32.70.Cs 



I. INTRODUCTION 

It is well known that Random Matrix Theory is able to 
reproduce successfully statistical properties which char- 
acterize the semiclassical regime of quantum systems 
with a few freedoms having a fully ergodic classical phase- 
space. This is true for bound systems since the short- 
ranged spectral statistics, those using energy levels which 
are separated by some multiples of the mean level spac- 
ing, are well described in this approach (see, e.g., Rcf. 
0). This is also true for scattering systems since the 
fluctuations of the S-matrix elements are well modeled 
with the help of random matrices 0). The predictions 
of Random Matrix Theory are universal, that is they 
are independent of the microscopic details of the partic- 
ular system under study. It is Classical Mechanics which 
provides the ultimate justification for universality. For 
instance, the universal behavior of the spectral statistics 
is justified by resorting to the principle that the very 
long unstable periodic orbits of the underlying classical 
dynamics are uniformly distributed in phase-space (see, 
e.g., Ref. ||). For instance also, the universal energy de- 
pendence of the ^-matrix autocorrelation function is ex- 
plained by the fact that the distribution in length of the 
very long unstable scattering trajectories has a generic 
functional form for classical ergodic systems (see, e.g., 
Ref. 1). 

Another example of universal behavior has to do with 
the transition matrix elements associated to an opera- 



tor perturbing a quantum system with an ergodic phase- 
space. Indeed, Random Matrix Theory predicts that 
these transition matrix elements are independent ran- 
dom variables distributed in a normal way ||. However, 
mean and variance characterizing these distributions are 
free parameters and therefore lack a physical meaning 
in this approach. In view of the previously given exam- 
ples, it appears natural to link this other universal prop- 
erty to Classical Mechanics again. This can actually be 
done by means of the generalization of the semiclassical 
theory pioneered by Gutzwiller to arbitrary (well 
behaved) operators K|, a generalization which allows 
to give a classical interpretation of mean and variance. 
Indeed, in the semiclassical limit of Quantum Mechan- 
ics, the leading order expression of the mean value of the 
diagonal transition matrix elements is equal to the micro- 
canonical average of the Weyl transform of the perturb- 
ing operator [pH. In the same limit, the leading order 
expression of the variance associated to the Gaussian dis- 
tributions generated with non-diagonal transition matrix 
elements is proportional to the Fourier transform of the 
classical autocorrelation function of the Weyl transform 
of the perturbing operator |^|| . The link between classi- 
cal variance and autocorrelation function has been stud- 
ied numerically in three quantum systems with an ergodic 
phase-space, a quartic oscillator [ |10[ , an oval billiard JO] 
and the Hydrogen atom in a strong magnetic field Jl2| . 
In these systems and for several examples of perturbing 
operators, it has actually been found that the autocorre- 
lation function of the Weyl transform of the perturbing 
operator determines with an excellent accuracy the lo- 
cal values of the classical variance. Moreover, it has been 
verified jll| that the semiclassical corrections to the lead- 
ing order expression of the variance which are due to the 
shortest unstable periodic orbits are in very good quanti- 
tative agreement with the theoretical predictions of Refs. 
H[| . The autocorrelation function has also been studied 
in the framework of a purely regular billiard for two ex- 
amples of perturbing operators |l3). It has to be added 
that a method has been proposed H] to compute the 
variance associated to the distribution of diagonal tran- 
sition matrix elements, a quantity which has not been 
considered in Refs. This method, which rests on pe- 

riodic orbit theory, has an interesting formal relationship 
with the supersymmetric methods used in the description 
of disordered systems . 

Generic quantum systems with a small number of free- 
doms d do have a mixed regular-ergodic phase-space. 
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This means that two different kinds of orbits appear in 
the underlying classical dynamics. There are orbits which 
wind regularly round d-dimensional tori and there are or- 
bits which explore densely (2d — l)-dimensional regions 
of the energy shell in a highly chaotic manner. Phase- 
space is therefore naturally partitioned into regular com- 
ponents, each including infinitely many neighboring tori, 
and ergodic components which are free of tori. These var- 
ious components are mutually independent. Since Ran- 
dom Matrix Theory appears to be relevant to quantum 
systems with a fully ergodic phase-space exclusively, the 
previously listed universal properties do not hold any- 
more for quantum systems with a mixed phase-space. 
However, one may still wonder about the presence in sta- 
tistical properties pertaining to mixed quantum systems 
of fingerprints related to the coexistence of regular and 
ergodic regions in phase-space. Berry and Robnik fl(f | 
have been among the first to tackle this question. Assum- 
ing that the semiclassical spectrum of a mixed quantum 
system is the superposition of statistically independent 
sequences of levels from each of the phase-space compo- 
nents, they have deduced a semi-phenomenological closed 
formula for the level spacing distribution. The formula 
makes use of the total phase-space volume of all regular 
regions as well as of the individual phase-space volumes 
of all ergodic regions. Their line of argument has been ex- 
tended to other spectral statistics |l7]] . The validity of the 
semi-phenomenological formulae obtained for the spec- 
tral statistics has been supported by numerical computa- 
tions |l7j . Apart from the one of Berry and Robnik jl6) , 
two other parametrizations of the level spacing distribu- 
tion have been proposed to characterize a mixed quan- 
tum system. These are the parametrizations of Brody 
fl8f and of Izrailev , whose validity has also been sup- 
ported by numerical computations. The link between the 
mixed character of phase-space and statistical properties 
of transition matrix elements has been studied by Rob- 
nik and Prosen |20pl| ]. Extending to eigenvectors the 
assumption of independent sequences originally proposed 
for energy levels, they have argued that the only transi- 
tion matrix elements which have to be considered in the 
semiclassical limit are those whose initial and final states 
can be associated to the same component of phase-space. 
Moreover, they have generalized the results of Refs. [^|J|] 
by showing that the expression of the variance character- 
izing the distributions of non-diagonal transition matrix 
elements whose initial and final states can be associated 
to a given component of phase-space (which is either an 
ergodic region or some set of quantized tori lying closely 
together) always involves a microcanonical average over 
this very component. Numerical computations [ p0|j2l| 
have corroborated the validity of their results. 

The purpose of the present paper is to provide a de- 
tailed study of the quantitative contributions coming 
from the various components of a mixed phase-space to 
the mean characterizing the distribution of diagonal tran- 
sition matrix elements and to the variance characterizing 
the distributions of non-diagonal transition matrix ele- 



ments. It means in particular to complement the study 
of Robnik and Prosen on a specific aspect of the link 
between the mixed character of phase-space and statis- 
tical properties of transition matrix elements which was 
not considered by them. This aspect concerns the use 
of the stable periodic orbits in the computation of the 
quantitative contributions coming from the various reg- 
ular components of phase-space to the values of mean 
and variance. The study of this particular point is in- 
teresting because of the two following reasons. On the 
one hand, the method proposed in Refs. [^[H) to com- 
pute the variance is not easy to use in practice since it 
requires the prior identification in phase-space of the ad- 
equate neighboring quantized tori. On the other hand, 
a given stable periodic orbit is in a way able all alone 
to account for the quantitative contribution of its asso- 
ciated regular component to the semiclassical density of 
states p^j2^ ,[7[|. Considering these two facts, one may 
wonder whether stable periodic orbits (whose practical 
identification is usually easier than the one of specific 
quantized tori) cannot help to compute in a convenient 
way the quantitative contributions of the regular regions 
to the values of mean and variance. The paper shows 
that this is indeed the case by providing closed formulae 
for mean and variance in which stable periodic orbits are 
involved. Such formulae generalize to the case of mixed 
quantum systems those studied in Refs. [p|Js| Jll| , [T2] | . This 
work extends the one of Ref. where details have been 
omitted because of lack of space. The illustrative numer- 
ical computations are again done in the framework of the 
Hydrogen atom in a strong magnetic field. 

The paper is organized as follows. In Section II, some 
details relating to the scaling property which character- 
izes the Hydrogen atom in a strong magnetic field are 
given and the choice of the perturbing operator used for 
the numerical computation of the transition matrix ele- 
ments is justified. Section III deals with the mean value 
of the distribution of diagonal transition matrix elements, 
a quantity which was not considered in Refs. |i0|]2l[ |. The 
quantitative contributions of the ergodic and regular re- 
gions making up the mixed phase-space are identified. 
The contributions of the various regular components are 
calculated with the help of the associated stable peri- 
odic orbits. The illustrative numerical computations are 
done for three different values of the scaled energy. Sec- 
tion IV is concerned with the variance characterizing the 
distributions of non-diagonal transition matrix elements. 
First, the orders of magnitude of the transition proba- 
bilities whose initial and final states are referring or not 
to the same phase-space component are compared. The 
quantitative contributions to the variance of the ergodic 
and regular regions are then studied. The manner in 
which a given stable periodic orbit is used in the cal- 
culation of the contribution coming from the associated 
regular component is discussed. The illustrative numer- 
ical computations are done for the same values of the 
scaled energy as in the previous section. A Conclusion 
summarizes the main results of the paper. 
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II. SCALING AND TRANSITION MATRIX 
ELEMENTS 

To begin with, it is not an easy task to study the quan- 
titative effect of the underlying classical dynamics on 
the statistical properties of transition matrix elements. 
This is due to the fact that the structure of phase-space 
changes with energy for most systems. Consequently, any 
diagonal transition matrix element (m|A|m) (|m) being 
an eigenvector of the system and A the perturbing oper- 
ator) is usually going with a different energy shell. The 
problem becomes even trickier when considering non- 
diagonal transition matrix elements (m|^4|n) since initial 
and final eigenstates are always going with different en- 
ergy shells. The study should therefore be done using a 
set of transition matrix elements in which initial and final 
eigenenergies are restricted to an energy interval small 
enough for phase-space to keep its structure practically 
unchanged. However, in order to get a significant number 
of transition matrix elements, one should choose an en- 
ergy region which is high enough in the spectrum so that 
the density of states can make up for the smallness of the 
used energy interval. This would obviously require the 
diagonalization of very large matrices in practice. Such 
numerical difficulties can be avoided by studying scaling 
systems. These are systems possessing scaling properties 
which imply that phase-space has the same structure at 
all energies. The Hydrogen atom in a strong magnetic 
field [£5| and the three-body Coulomb system |26) are 
physically interesting examples of scaling systems. The 
Hamiltonian of a scaling system depends on a scaling 
parameter. It is useful to follow the variation of the en- 
ergy levels with the scaling parameter p5p6| ] since one 
can extract informations about the underlying classical 
dynamics from spectra taken at different values of the 
scaling parameter by application of the so-called method 
of scaled energy spectroscopy |^,^8| . 

In atomic units, the quantum Hamiltonian of the Hy- 
drogen atom in a strong magnetic field reads 

H 7 (p, f ) = ip 2 - I + \ 1 L Z + i 7 2 (£ 2 + f) . (1) 

Here L z is the component of the angular momentum op- 
erator along the direction of the magnetic field. This 
component is conserved and, consequently, the azimuthal 
quantum number m is a good quantum number. Numer- 
ical computations are restricted to the subspace m = 
in this study. The Hamiltonian is invariant under reflec- 
tion with respect to the plane which is perpendicular to 
the direction of the magnetic field and the z-parity n z 
is thereby a good quantum number too. The param- 
eter 7 is the magnetic field strength in atomic units, 
7 = £>/(2.35 x 10 5 T). When expressed in terms of 
the scaled coordinates f = 7 2 / 3 r and scaled momenta 
p = 7~ 1 / 3 p, the classical Hamiltonian scales as 

ff 7 =i(p, f ) = ip 2 - ~ + \l z + \{i 2 + f) 



7 



- 2 / 3 # 7 (p,r) = 1 ~ 2 / 3 E = E 



(2) 



with E the excitation energy Therefore, the classical 
dynamics obtained from the scaled equations of motion 
does not depend on E and 7 independently but on a 
single parameter combining both physical quantities, the 
scaled energy E — j~ 2 / 3 E. This implies that the struc- 
ture of phase-space is identical for any pair (E,j) leading 
to the same value of the scaled energy. 

If z and p z are the scaled coordinate and momentum 
operators along the direction of the magnetic field, one 
has [z,p z ] = i^^Ti by virtue of the previously given 
definition of scaled variables. The dependence of the 
quantum dynamics on the magnetic field strength 7 can 
thereby be taken into account by means of an effective 
Planck's constant h c s = 7 1 ^ 3 ?i (h = 1 subsequently). 
One approaches the semiclassical limit at constant scaled 
energy ?i ff — > by decreasing the value of 7. The quan- 
tization of Eq. (g) in the subspace m = leads to a 
generalized eigenvalue equation for the scaling parame- 
ter w — 7 -1 / 3 = h~g. The introduction of the scaled 
semi-parabolic coordinates p = \Jf + z and v = y/f — z 
allows to write this generalized eigenvalue equation as 



2E{p 2 + D 2 )- -^V(m 2 + *> 2 )+4 



w 



Pu+Pu) ^(A,^) 



(3) 
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with the radial operators p^ and p v defined as 
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v dv 



dv 



Eq. (^) can be written in matrix form by using the com- 
plete set of basis functions which is composed of the ten- 
sorial products of the eigenvectors of two uncoupled two- 
dimensional harmonic oscillators with frequency y—2E 
p5| , p9| . The generalized eigenvalue equation is solved 
with the help of the Lanczos spectral transformation 
method, which is adapted to the diagonalization of sparse 
symmetric matrices |50|. One obtains in this way the 
spectrum of eigenvectors |\& m ) and corresponding eigen- 
values w m = ( 7 _1 / 3 ) m . The structure of phase-space is 
the same at every eigenenergy w m since the generalized 
eigenvalue equation has been solved at constant scaled 
energy E. Consequently, the matrix elements describing 
the various transitions caused in the spectrum by a per- 
turbing operator are all going with the same underlying 
classical dynamics. It has to be noted that the eigen- 
vectors I "I'm) are not orthogonal. However, the modified 

/),, +p~\^ m ) (with the same corre- 
are orthogonal, i.e. 



eigenvectors \m) - 
sponding eigenvalues w r , 



(* m |p!+2Vl*n) =S„ 



(4) 



The perturbing operator which has been chosen for this 
study is 
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A = 



w 



2rp 



~2 



x 2 



(5) 



The second expression in this equation is the effective ex- 
pression of the operator A when acting onto eigenvectors 
belonging to the subspace m — 0. The same operator 
has been used in previous papers [JTJ 24 already. The 
reason is that this particular operator is convenient to 
handle in practice since the general expression of the as- 
sociated transition matrix elements in the orthonormal 
basis {|m)} of modified eigenvectors is 



(m\A\n) = w m w n (^ m \^n) 



(6) 



The computation of the transition matrix elements at 
constant scaled energy amounts therefore merely to the 
computation of the various overlaps of the eigenvectors 
{l^m)}- Moreover, this computation can be restricted to 
a subspace which is labeled by a given value of tt z because 
the operator A connects only eigenstates of the same z- 
parity. Numerical computations are done here in the sub- 
space m* = + . This special perturbing operator offers 
another computational advantage. Indeed, it turns out 
that the average at constant^ scaled energy of the Weyl 
transform A of the operator A along a given trajector y in 
phase-space has an analytic expression, given in Eq. jlOj ) 
below. This average, which is important for the phys- 
ical interpretation of the numerical results presented in 
this paper, can therefore be computed accurately. The 
Weyl transform of the perturbing operator is expressed 
most simply with the help of the scaled semi-parabolic 
momenta p^ and p v , the conjugate variables of the scaled 
semi-parabolic coordinates p and v. These momenta are 
defined as the derivatives of the coordinates jl and v with 
respect to the so-called rescaled time variable r ]2j| . The 
expression of A is 



1 



Pi 



■Pi 



(7) 



The scaling parameter w = and the scaled ac- 

tion s(E) = J (p^djt + p v dv) are conjugate variables 
at constant scaled energy. Since the variable w is an 
energy variable, the variable s can be considered as a 
time variable which measures the length of trajectories in 
the four-dimensional phase-space spanned by the scaled 
semi-parabolic variables (jl, v-,Pix,pv 

) H). The infinites- 
imal scaled action ds and the infinitesimal rescaled time 
dr are linked together by the relation 



ds = p^djl + p v di> 



(Pi 



Pv)dr 



(8) 



This relation implies that the Weyl transform A can be 
expressed at constant scaled energy as the derivative of 
the variable r with respect to the variable s, i.e. 



Aip^pv) 



dr 
ds 



A(s) 



(9) 



The average at scaled energy E of the Weyl transform A 
along a trajectory with scaled action S(E) corresponding 
to a total rescaled time r is then simply 



1 



S(E) 



III. 



S(E) 



dsA(s) = 



1 



S(E) 



dr = 



S(E) 



(10) 



DIAGONAL TRANSITION MATRIX 
ELEMENTS 



This section studies the quantitative contributions 
coming from the ergodic and regular components of the 
mixed phase-space to the mean value characterizing the 
distribution of diagonal transition matrix elements asso- 
ciated to the perturbing operator A. The study is done 
for three successive values of the scaling energy E, which 
are E = -0.2, E = -0.316 and E = -0.4. Fig. | displays 
the Poincare surfaces of section for these different cases 
in the scaled semiparabolic representation v = 0). 

In the case E = —0.2 (Fig. [TJa), the surface of section has 
two islands of stable motion which are embedded in a 
simply connected ergodic sea. Both islands belong to the 
same regular component since they are conjugate with 




FIG. 1. Poincare surfaces of section of the Hydrogen atom 
in a strong magnetic field in the scaled semi-parabolic repre- 
sentation {fi,Pn \ v — 0) for three different values of the scaled 
energy E. (a) E = -0.2; (b) E = -0.316; (c) E = -0.4. 
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FIG. 2. Distribution of values of the diagonal transi- 
tion matrix elements (m|yl|m) corresponding to the operator 
A — l/2rp 2 in the case of the scaled energy E — — 0.2. (a) 
Distribution of values of the matrix elements associated to 
the chaotic (x) and regular (+) states as a function of the 
scaling parameter w; (b) Distribution of values of the matrix 
elements associated to the regular states as a function of the 
inverse of the mean total density of states po,t(w). 

respect to the line p^ — 0. This regular component is 
associated to a stable periodic orbit lying in the plane 
perpendicular to the direction of the magnetic field, with 
scaled action S r ~ 6.49086. Fig. ||a displays, for the 
same value of the scaling energy, the distribution of val- 
ues of the diagonal transition matrix elements as a func- 



This 



tion of the scaling parameter w = j' 1 -' 6 = fi~J . 
figure shows clearly the presence of two different pat- 
terns. On the one hand, there appears a statistical dis- 
tribution built up with diagonal transition matrix ele- 
ments which are represented by "x" symbols. A further 
analysis reveals that these matrix elements involve eigen- 
states whose Husimi distribution ||^] is essentially local- 
ized within the ergodic sea. Such eigenstates are called 
"chaotic" in the sequel. They are generically labelled as 
|m c ), with corresponding eigenvalues w$ . These eigen- 
values belong to a subset of the spectrum which is charac- 
terized by a mean density of states po iC (w) — 0.702u>. On 
the other hand, one distinguishes in the lower part of Fig. 
||a several sequences of diagonal matrix elements which 
are represented by "+" symbols. The Husimi distribu- 
tions of the corresponding eigenstates are well localized 



within the two islands of stable motion in Fig. [j]a. Such 
eigenstates are called "regular" in the sequel. They are 
generically labelled as \m r ), with corresponding eigen- 

(r) 

values Wm ■ These eigenvalues belong to a subset of the 
spectrum which is complementary to the previous one 
and whose mean density of states is po. r (w) = 0.066w. 
The coexistence of these two different patterns of diag- 
onal transition matrix elements illustrates well the fact 
that the spectrum of a mixed quantum system can be 
split into independent subsets of eigenvectors, each of 
which going with a different component of phase-space 

& • ■ B - - 

The quantitative interpretation of Fig. G& is done with 

the help of a scaled spectral function D(w) which takes 
all diagonal matrix elements of the perturbing operator 
A into account. It is defined as 



D{w) 



(m\A\m)S(w — iv m ) 



(11) 



According to Refs. f||§, the semiclassical expression of 
D(w) which is relevant to the fully ergodic case is 



D(w) 



dpndp v dp,du (jl 2 + i> 2 



x A{p ll ,p v )8{H{p ll ,p lJ ,jl,i)) - 2) 



00 sin s 



[det(M* - 1)] 



1/2 



(12) 



with 



1 

2 



pl+Pl)-(t 2 + i> 2 )E+^ 2 v 2 tf + v 2 ). (13) 



The expression of D(w) is written here in the scaled semi- 
parabolic representation of phase-space. The first term 
on the r.h.s. of Eq. (|l2|), the so-called Weyl term, gives 
the classical contribution to the scaled spectral function. 
It factorizes into two parts as a consequence of the use of 
the scaled semiparabolic representation. One part de- 
pends solely on the scaling parameter w whereas the 
other part, which involves a phase-space integration over 
the whole energy shell O at constant scaled energy E, 
is independent of it. The second term gives the leading 
order corrections to the Weyl term in an asymptotic ex- 
pansion of D(w) into powers of Ti. These corrections are 
generated by the periodic orbits of the underlying clas- 
sical dynamics, which are all unstable. The outer sum 
in the second term runs over the set {p} of all primitive 
periodic orbits whereas the inner sum runs over all rep- 
etitions s of every primitive periodic orbit. Each primi- 
tive periodic orbit p is characterized by an amplitude A p 
whose expression is 



A„ 



dsA(s) 



(14) 
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In this equation, the integration is done along the con- 
sidered periodic orbit and over one period S P (E) = 
§(p^dp + p v dv). The scaled action S P (E) appears in 
the argument of the sine function together with the scal- 
ing parameter w since it scales as S p (E) = w~ 1 S p (E), 
S P (E) being the usual action computed at the excita- 
tion energy E = E/w 2 . Every primitive periodic orbit is 
also characterized by its monodromy matrix M p and its 
Maslov index a p , which are both independent of energy 
for scaling systems. The importance of the correction 
due to the s-th repetition of the primitive periodic orbit 
p depends on the determinant of the difference between 
the s-th power of M p and the unit matrix /. As Fig. |2|a 
shows, the scaled spectral function D(w) can be written 
as the sum of two different contributions in the case of a 
mixed system, i.e., 



D(w) = D c (w)+D r (w) 



(15) 



with 



%) = J2 K(r)|i|m c(r) ) 5{w - wt [r)) ) ■ (16) 



The semiclassical expression of D c (w) (resp. D r (w)), the 
scaled spectral function which is associated to the ergodic 
(resp. regular) component of phase-space, is similar to the 
one given in Eq. (|l^) . The integration in the Weyl term is 
now done over the corresponding part of the energy shell 
at constant scaled energy. The leading order corrections 
to the Weyl term are due to the unstable (D c {w)) or 
stable (D r (w)) periodic orbits of the underlying classical 
dynamics. 

The semiclassical expression of D c (w) predicts that the 
classical contribution to the mean value (A) c which char- 
acterizes the distribution of diagonal transition matrix el- 
ements corresponding to the chaotic eigenstates is given 
by the microcanonicai average of the operator A over the 
ergodic part fl c of the energy shell i.e., 



J n dp^dp v dpdv (p 2 + v 2 ) A{p,,,,p v )8{H - 2) 



J n dp^dpudjldv (jl 2 + v 2 ) 5(H 



(17) 



with H(p^p u ,p,D) defined in Eq. (|3|) |3|. This for- 
mula is not easy to use in practice for the computa- 
tion of the classical contribution since it requires the nu- 
merical calculation of three-dimensional phase-space in- 
tegrals. However, just as in Refs. the ergodic the- 
orem allows to replace the microcanonicai average with 
a time (i.e. scaled action) average of the Weyl transform 
of the perturbing operator along any trajectory explor- 
ing the ergodic part of the energy shell in a uniform way. 
With the help of Eq. (|l0|), the classical contribution to 
the mean value (A) c can be finally expressed as 



lim ^ / dsA(s) 



t(S) 



S^oo S 



t(S) being the total rescaled time corresponding to the 
scaled action S. The ergodic theorem ensures that the 
limit on the r.h.s. of this equation is well defined and 
unique. This formula is well suited for the numerical cal- 
culation of the classical contribution to the mean value. 
In the case E — —0.2, it gives the value (A) c — 0.41, 
which is in good agreement with the value 0.411 obtained 
from the statistical distribution of Fig. ||a. Eq. ( |l8|) can 
therefore be used to compute with an excellent precision 
the mean value characterizing the distribution of diago- 
nal matrix elements associated to the chaotic eigenstates. 

Fig. |]a suggests that all sequences of diagonal matrix 
which correspond to the regular eigenstates are converg- 
ing towards the same limit as the scaling parameter w in- 
creases. This is indeed the case, as clearly shown by Fig. 
||b. In this figure, the values of the diagonal matrix ele- 
ments pertaining to the seven sequences of the previous 
figure are represented as a function of the inverse of the 
mean total density of states po,t(w) = pa. c (w) + po. r (w). 
Such a representation is useful since it allows to extrapo- 
late the values of the individual diagonal matrix elements 
in the semiclassical regime. This comes from the fact 
that the mean total density of states is proportional to 
w in the case of the Hydrogen atom in a strong magnetic 
field (po,t = 0.768u; in the case E = —0.2), and so its 
inverse is proportional to the effective Planck's constant 
freff = 7 1 ^ 3 of this system. The semiclassical limit at 
constant scaled energy h e s — » is therefore reached as 
l/po,t(w) — > 0. One sees that the values of the diago- 
nal matrix elements in the sequences are approximated 
all the better by a common limiting value (A) r as one 
comes closer to the semiclassical regime. This behavior 
is in sharp contrast with the one of the diagonal matrix 
elements corresponding to the chaotic eigenstates, which 
are distributed statistically around the mean value (A) c . 
The limiting value (A) r can be most easily computed by 
using an expression of the scaled spectral function D r (w) 
which is obtained by resumming formally the leading or- 
der corrections to the Weyl term generated by the stable 
periodic orbit and all its repetitions in the following way 
p3|,p|J^] . The expression of these corrections for a partic- 
ular stable periodic orbit r is ppM^ 



i y—^ exp(is(wS r (E) — a r n/2)) 

— Ti r y ; r 

7T ^— ' 2sin(s7T7 r ) 



(19) 



Here j r is the winding number of the stable orbit, which 
is independent of energy for scaling systems. After ex- 
pansion of the sine function in the denominator and the 
subsequent use of the Poisson summation formula, this 
expression can be rewritten as 



A r J2 5 ( wS r (E) -2tt 



■In • Ar — — j 



(20) 



lim 



(18) 



This is a formal expression of the scaled spectral func- 
tion D r (w) in the semiclassical regime. It predicts that 



G 



the semiclassical spectrum of the regular states is similar 
to the spectrum of a two-dimensional harmonic oscilla- 
tor 1 23 g,0]. One direction of harmonic motion is along 



the stable periodic orbit, whereas the other direction of 
harmonic motion is transverse to it. This motion takes 
place on neighboring tori surrounding the orbit. Each 
such torus is identified by two quantized scaled actions, 
the action /„ = 27r(n + a r /4) which is associated to the 
longitudinal motion and the action Ik = 2~K(k + l/2)7 r 
which is associated to the transverse motion [pi| . Ev- 
ery regular state is therefore labelled by two quantum 
numbers n and k. The range of values of both quan- 
tum numbers is bounded in practice as a consequence of 
the finite size of the corresponding region of stable mo- 
tion J35| . It has been checked that all regular eigenstates 
building up a given sequence of diagonal matrix elements 
in Figs. [|a,b are labelled by the same quantum number 
k. This quantum number is even as a result of the par- 
ity symmetry with respect to the plane perpendicular to 
the direction of the magnetic field. The lowest sequence 
corresponds to k = 0, the next one to k = 2 and so on 
until the last identified sequence which corresponds to 
k = 12. Fig. ||b shows that there are a few diagonal ma- 
trix elements which do not belong to any sequence. It 
has been checked that all of them correspond to regular 
eigenstates which are involved in quasi-crossings and so 
cannot be labelled by a fixed value of k. The quantum 
number n differentiates the regular eigenstates which are 
associated to a given sequence. The previous formal ex- 
pression of the scaled spectral function D r (w) predicts 
also that the limiting value (A) r of the sequences is pro- 
portional to the amplitude A r of the stable periodic orbit, 
whose expression is given by Eq. (|l4| ) . The amplitude A r 
has to be divided by the scaled action S r of the orbit for 
the purpose of normalization. With the help of Eq. (10), 
this limiting value can be expressed as 



(A) r = ^- U~sA(s) = £ 



(21) 



T r being the total rescaled time corresponding to the 
scaled action S r . It has to be remarked that the classi- 
cal contribution (A) c to the mean characterizing the dis- 
tribution of diagonal matrix elements associated to the 
chaotic eigenstates, Eq. (|l8|) , and the classical limiting 
value (A) r of the sequences of diagonal matrix elements 
corresponding to the regular eigenstates, Eq. (|2l|), have 
similar expressions. Each expression makes use of the 
classical trajectory which is related to the component of 
phase-space going with the studied subset of diagonal 
matrix elements. In the case E = —0.2, the scaled action 
corresponds to a total rescaled time r r — 2.4135. This 
gives the value (A) r — 0.372 which, as seen in Fig. ||b, 
agrees very well with the value obtained as the common 
intersection of the dashed lines interpolating the values 
of the diagonal matrix elements as one goes towards the 
semiclassical limit. 

Since the total distribution of diagonal matrix elements 



is composed of two different subsets corresponding to the 
two different components of phase-space, its mean value 
(A} t is expected to be the sum of two properly weighted 
contributions. The contribution which comes from the 
mean value (A) c (resp. limiting value (A) r ) is weighted 
by the ratio of the mean density of states of the subset 
of chaotic (resp. regular) states to the mean total density 
of states. Each weighting factor is roughly equal to the 
relative volume of the corresponding component of phase- 
space. The expected expression of (A) t is therefore the 
following 



\Po,t{ w )J \Po,t(w)J 



(22) 



In the case E = —0.2, this formula gives the value (A) t — 
0.407, which is very close from the value (A) t = 0.409 
obtained from the whole distribution of Fig. ||a. Con- 
sequently, Eq. ( p2[ ) allows to calculate with an excellent 
precision the mean value characterizing the total distri- 
bution of diagonal matrix elements. It generalizes to the 
case of a scaled system with a mixed phase-space the ex- 
pression of the mean value given in Refs. ||D for a scaled 
system with an ergodic phase-space. 

It is expected that there are as many limiting values 
(A) r as there are regular components in phase-space. 
This important point can be checked by studying the 
distribution of diagonal matrix elements for values of the 
scaled energy for which several stable periodic orbits are 
present in the underlying classical dynamics. This is done 
in the cases E = —0.316 and E = —0.4. The Poincare 
surface of section corresponding to the first case is dis- 
played in Fig. |l|b. There are five islands of stable mo- 
tion on both sides of the line p M = 0, a central island 
surrounded by a chain of four islands. The two conju- 
gate central islands belong to the same regular compo- 
nent which, as in the case E = —0.2, is associated to 
a stable periodic orbit lying in the plane perpendicular 
to the direction of the magnetic field. This orbit has a 
scaled action S r i — 7.903534, corresponding to a total 
rescaled time r r \ — 3.108460. The shape of this orbit 
is depicted in the semi-parabolic coordinate representa- 
tion (p, i>) in the lower inset of Fig. ||b. The two conju- 
gate chains of four islands of stable motion surrounding 
the central islands belong to another regular component. 
The associated stable periodic orbit has a scaled action 
S r 2 = 12.158133, corresponding to a total rescaled time 
r r 2 = 5.028604. The shape of this orbit is depicted in the 
scaled semi-parabolic coordinate representation (/2, u) in 
the upper inset of Fig. ||b. Fig. ||c shows the Poincare 
surface of section corresponding to the case E = —0.4. 
As in the two previous cases, the two conjugate islands of 
stable motion on both sides of the line = belong to 
a regular component with a stable periodic orbit lying in 
the plane perpendicular to the direction of the magnetic 
field. This orbit has a scaled action S r i = 7.02481 corre- 
sponding to a total rescaled time T r \ — 2.866. The stable 
island in the middle of the surface of section pertains to 
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FIG. 3. Same as Fig. 2 but in the case of the scaled 
energy E — —0.316. The shapes of the two stable periodic 
orbits which are used in the analysis of this case are shown in 
the insets. 

another regular component, with a stable periodic orbit 
which is parallel to the direction of the magnetic field. 
This orbit has a scaled action S r 2 = 5.791216, corre- 
sponding to a total rescaled time r r 2 = 2.896. 

Fig. ||a (resp. Fig. |Ja) displays the distribution of val- 
ues of the diagonal matrix elements as a function of 
the scaling parameter w in the case E = —0.316 (resp. 
E = —0.4). Just as in the case E = —0.2, the whole 
distribution results from the juxtaposition of a statisti- 
cal distribution connected to the ergodic component of 
phase-space and a set of sequences connected to the reg- 
ular components of phase-space. These sequences fall 
into two different groups, as the fingerprint of the exis- 
tence of two different stable periodic orbits. This is par- 
ticularly clear in Fig. ^a, where each group of sequences 
is on a different side of the statistical distribution. In 
the case E = —0.4, the long streaks of diagonal matrix 
elements appearing in the statistical distribution are pre- 
sumably the mark of one or several other regular compo- 
nents, which are associated to the small islands of stable 
motion surrounding the three main islands in Fig. |l|c. 
As previously explained, the identification of the group 
of sequences linked to a particular regular component is 
most easily done by representing the values of the diago- 
nal matrix elements in all sequences as a function of the 
inverse of the mean total density of states poj(w). This 
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FIG. 4. Same as Fig. 2 but in the case of the scaled energy 
E = -0.4. 

is done in Fig. ||b (resp. Fig. ib) for the case E — —0.316 
(resp. E = —0.4). In Fig. 3 b, one sees that the three 
lower and the five upper dashed lines interpolating the 
values of the diagonal matrix elements as one goes to- 
wards the semiclassical limit have a common intersection. 
The value of the common intersection of the three lower 
sequences is very precisely equal to the limiting value 
(A) r i — T r \/S r \ — 0.3933 of the stable periodic orbit 
lying in the plane perpendicular to the direction of the 
magnetic field. The regular component which is related 
to these sequences is therefore identified in a straight- 
forward way. It is obviously the one to which the two 
conjugate central islands of stable motion of Fig. |l|b are 
pertaining. As in the case E = —0.2, the three lower se- 
quences are labelled by even values of the quantum num- 
ber k, with the lowest sequence corresponding to k = 0. 
The value of the common intersection of the five upper 
sequences is very precisely equal to the limiting value 
{A) r2 — T r i/S r 2 — 0.4136 of the other stable periodic 
orbit. Consequently, the regular component which is as- 
sociated to these sequences is the one whose fingerprint 
in Fig. |l|b is constituted by the two conjugate chains of 
four islands of stable motion. Since this other orbit has 
no symmetry with respect to the plane perpendicular to 
the direction of the magnetic field, the five upper se- 
quences are labelled by both even and odd values of the 
quantum number k, with the uppermost sequence cor- 
responding to k = 0. Similar results are obtained in 
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the case E = —0.4. Indeed, one sees in Fig. |b that 
the groups of dashed interpolating lines corresponding to 
the sequences above and below the statistical distribu- 
tion have each a common intersection. The value of the 
common intersection of the sequences below the statis- 
tical distribution is very precisely equal to the limiting 
value (A) r i — T r \/S r i = 0.408 of the stable periodic or- 
bit lying in the plane perpendicular to the direction of 
the magnetic field. The regular component connected to 
these sequences is, as in the two previous cases, the one 
to whom the two conjugate islands of stable motion in 
the Poincare surface of section are belonging. This means 
that the limiting value of a group of sequences readjusts 
itself to the change in length of the associated orbit as 
one moves from a particular example of the mixed phase- 
space to another. The sequences below the statistical dis- 
tribution are labelled by even values of k 7 with the lowest 
sequence corresponding to k = 0. As expected, the value 
of the common intersection of the sequences above the 
statistical distribution is (nearly) equal to the limiting 
value (A) r 2 = T r ijS r i = 0.500 of the stable periodic or- 
bit which is parallel to the direction of the magnetic field. 
By way of consequence, the regular component which is 
linked to these sequences is the one whose fingerprint in 
Fig. 0c is the central island of stable motion. These se- 
quences are labelled by both even and odd values of k, 
with the uppermost sequence corresponding to k = 0. 
Finally the formula in Eq. (|2^) giving the mean value 
which characterizes the complete distribution of diago- 
nal matrix elements is generalized to the case of several 
regular components simply by adding as many weighted 
classical values (A) r as there are stable periodic orbits in 
the underlying classical dynamics. 



IV. NON-DIAGONAL TRANSITION MATRIX 
ELEMENTS 

This section studies the quantitative contributions of 
the ergodic and regular components of the mixed phase- 
space to the variance characterizing the distributions of 
non-diagonal transition matrix elements associated to the 
perturbing operator A. This is done with the help of 
scaled spectral functions which take the different types 
of non-diagonal transition matrix elements into account, 
those coupling chaotic or regular eigenstates together as 
well as those coupling a chaotic and a regular eigenstate 
together. Since they deal with non-diagonal matrix ele- 
ments, these scaled spectral functions depend necessarily 
on two energy scales, the scaled energy w and the scaled 
energy difference Aw. The scaled total spectral function 
C(w, Aw) is defined to be the sum of these scaled spectral 
functions, i.e., 



C(w, Aw) = C cc (w, Aw) + C rr (w, Aw) 
+ C cr (w, Aw) + C rc (w, Aw) 



(23) 



C a/ 3(w,Aw) 
= ^2 \(m a \A - (A) a S a p\np)\ 2 5, l 



Win +W n 



(24) 



and a = c, r; (3 — c,r. The (Lorentzian) smoothings 
of the ^-functions, of widths r\ and e, are introduced so 
that one can compute these spectral functions in prac- 
tice in spite of the discreteness of the spectrum. The 
values of the widths which have been chosen for the 
numerical computations are r/ = 5.0 and e = 0.02, in 
units of the scaling parameter w. The classical values 
(A) c , {A) r are subtracted from the appropriate spectral 
functions in order to eliminate the quantitative contribu- 
tion of the diagonal matrix elements in the semiclassical 
regime. The local variance cr^iw, Aw) which is associated 
to the statistical distributions of non-diagonal matrix el- 
ements coupling chaotic eigenstates together is related to 
the spectral function C cc {w, Aw) and to the mean den- 
sity of states po, c (w) by the formula B 



a^.(w, Aw) 



C cc (w, Aw) 

{poA w )) 2 



(25) 



By analogy with this formula, the local variance 
af(w,Aw) which is associated to the distributions of all 
non-diagonal matrix elements is defined through the ex- 
pression 



of (w, Aw) 



C(w,Aw) 
(po,t( w )) 2 



(26) 



with 



The quantitative contribution of each spectral function 
in Eq. ( p3| ) to the variance (T^(w,Aw) is studied in the 
sequel. 

The spectral function C cc (w, Aw) (resp. C rr (w, Aw)) 
is connected to the regular (resp. ergodic) component 
of phase-space in a manner which is detailed below. 
On the contrary, the spectral functions C cr (w, Aw) and 
C rc (w, Aw) cannot be associated to a particular compo- 
nent of phase-space since they take both chaotic and reg- 
ular eigenstates into account. According to Refs. pfj|j2l| ], 
the quantitative contributions of these two last spectral 
functions to C(w,Aw) should therefore be significantly 
smaller than those of the two first spectral functions. Fig. 
|5| shows that this is indeed the case in practice. This 
figure displays the distribution of values of the transi- 
tion probabilities |(n|A|m)| 2 from a given initial state 
|n) to a subset |m) of final eigenstates corresponding to 
a finite range of the energy spectrum at scaled energy 
E = — 0.2. In Fig. ||a, the chosen eigenstate \n) (the 
575 th state above the groundstate) is a regular eigenstate 

(r) 

with k = and eigenvalue w n ~ 38.5. On the contrary, 
the eigenstate \n) in Fig. (the 944 th state above the 
groundstate) is a chaotic eigenstate. In both figures, the 
diamonds (resp. crosses) are recording the values of the 
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FIG. 5. Distribution of values of the transition probabili- 
ties | (n | A | m) | 2 as a function of the scaling parameter w in the 
case of the scaled energy E = —0.2. (a) Regular eigenstate 
|n); (b) Chaotic eigenstate |n). The diamonds and crosses 
mark the transitions to regular and chaotic states \m), re- 
spectively. 

transition probabilities when the final eigenstate |m) is 
a regular (resp. chaotic) one. Each of the three visible 
sequences of diamonds in Fig. |5^l is associated to regular 
eigenstates which are labelled by the same (even) value of 
k, with the uppermost sequence corresponding to k — 0. 
This figure shows clearly that the probabilities of most 
transitions coupling regular eigenstates together are sev- 
eral orders of magnitude larger than those coupling reg- 
ular and chaotic eigenstates together. This implies in re- 
turn that the values of the spectral function C rr (w, Aw) 
are much larger that those of the spectral functions 
C rc (w, Aw) and C cr (w, Aw). Two other observations in 
Fig. [5|a are of interest. The first is that the square of 
the diagonal matrix element (n — 575|A|n = 575) is at 
least one order of magnitude larger than the probabili- 
ties of transitions coupling different regular eigenstates 
together. The second observation is that the probabili- 
ties of the transitions involving regular eigenstates with 
the same value of k are more than one order of magni- 
tude larger than those involving regular eigenstates with 
different values of k. Moreover, the larger the differ- 
ence in k between two regular eigenstates, the smaller 
the corresponding transition probability [M . It has also 
been checked that the values of the transition probabil- 
ities belonging to the sequence of diamonds with k = 



are decaying exponentially as one goes away from the 
diagonal matrix element. This type of decay is a univer- 
sal feature of one-dimensional systems Q . The system 
behaves therefore as an effective one-dimensional system 
at the level of the transition amplitudes which are con- 
nected to the first quantized torus surrounding the sta- 
ble orbit. Fig. ||b shows clearly that the probabilities 
of most transitions involving chaotic eigenstates are also 
several orders of magnitude larger than those involving 
chaotic and regular eigenstates. As a consequence, the 
values of the spectral function C' cc (w, Aw) are also much 
larger than those of the spectral functions C rc (w, Aw) 
and C cr (w, Aw). The quantitative contribution of the 
two last spectral functions to the variance af(w,Aw) is 
therefore negligible. 

As shown in Refs. ||||, the leading order contribution 
to the spectral function C cc (w, Aw) is proportional to the 
Fourier transform of the classical autocorrelation func- 
tion of the Weyl transform of the perturbing operator. 
By virtue of the ergodic theorem, the microcanonical av- 
erage over the ergodic part of the energy shell appearing 
in the leading order contribution can be replaced with a 
time (i.e. scaled action) average along any ergodic trajec- 
tory exploring this part in a uniform way. The resulting 
expression of the leading order contribution to the local 
^(w, Ai 



variance at. 



\w) is 



\(w, Aw) 



1 

irpo, c (w)' 



Re / dse i(Aw+le) 



C A, 



(5). 



(27) 

Here C^ c (s) is the classical autocorrelation function of 

the Weyl transform A(s), as computed along an ergodic 
trajectory of arbitrary large scaled action S, i.e., 



x 1(5' + 5/2) - (i) c U(5' - 5/2) - (A) 



(28) 



The ergodic theorem ensures that the limit on the r.h.s. of 
this equation is well defined and unique. The appearance 
in Eq. ( p7j ) of a damping factor containing the width 
e is a consequence of the smoothing of the ^-function 
which is associated to differences in scaled energies in 
Eq. (HJ). The formula in Eq. @) is well suited for the 
numerical computation of the classical contribution to 
the local variance a%(w, Aw). It shows that, at the level 
of the leading order contribution, the rescaled variance 
<Tc(w, Aw)po iC (w) is a function which depends only on the 
scaled energy difference Aw and no more on the scaled 
energy w. 

As a first numerical study, it is interesting to compare 
the values of the rescaled variances of (w, Aw)po,t(w) 
and o"e (if, Aw)po :C (w) over a large range of values of the 
scaled energy difference. This is done in Fig. ^ for the 
case E = —0.2. The full (resp. dashed) curve corresponds 
to the first (resp. second) rescaled variance. One sees that 
both curves are modulated by the period S r of the stable 
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FIG. 6. Rescaled variances of the distributions of 
non-diagonal transition matrix elements as a function of the 
scaled energy difference Aw in the case of the scaled en- 
ergy E = —0.2. Solid line: erf (w, Aw)/9o,t(w); Dashed line: 
<t;?(w, Aw)p 0lC (w). 

orbit lying in the plane perpendicular to the direction 
of the magnetic field since they exhibit peaks exactly at 
the integer multiples of the value Aw = 2-k/S t — 0.968. 
This modulation effect of the local variance was already 
pointed out in Rcfs. JI^JTl] ]. One observes also that the 
two curves are only slightly different from each other. 
This observation simply reflects the fact that the num- 
ber of transitions involving chaotic eigenstates is very 
much larger than the number of transitions involving reg- 
ular eigenstates in the chosen example. Consequently, a 
numerical analysis using the bare variances themselves 
would not allow to get precise enough values for the 
quantitative contribution to the total variance coming 
from the regular part of phase-space. As seen below, it 
is necessary to use the Fourier transform of the rescaled 
variances in order to extract quantitative precise results 
from the numerical data. In practice, the full curve in 
Fig. ^ has been calculated with the help of the numerical 
data built up from the distributions of the non-diagonal 
matrix elements whereas the dashed curve has been com- 
puted by using Eqs. It is therefore interesting to 
compare the predictions of the expression of the lead- 
ing order contribution to the local variance <j^.(w,Aw) 
with the numerical values obtained from the distribu- 
tions of non-diagonal matrix elements coupling chaotic 
eigenstates together. This is done in Fig. ^a for the case 
E = -0.2 and in Fig. |a for the case E = -0.316. In 
both figures, the full curve represents the Fourier trans- 
formC(S') of the rescaled variance a^(w, Aw)po tC (w), i.e., 



C(S) 



d(Aw) cos(AwS) <t c (w, Aw) po iC (w) , 

(29) 



as calculated from the relevant numerical data. On the 
other hand, the dashed curve represents the classical au- 
tocorrelation function C jj C {S), Eq. (|28|), as computed 
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FIG. 7. Fourier transform of the rescaled variances 
of the distributions of non- diagonal transition matrix ele- 
ments as a function of the scaled action S/2tt in the case 
of the scaled energy E — —0.2. (a) Solid line: Fourier 
transform of a^(w, Aw) po, c (v}); Dashed line: Autocorrela- 
tion function C \ c . (b) Solid line: Fourier transform of 
of (w, Aw)po,t{w); Dashed line: Autocorrelation function 
Cj c . (c) Solid line: Fourier transform of of(lu, Aw)po,t(w); 
Dashed line: Weighted sum of autocorrelation functions C ^ c 
and C4 r . 

along an arbitrary ergodic trajectory. Eq. ( |27| ) predicts 
that both curves are identical, i.e. C(S) = C ^ C (S) . Figs. 
0a and ||a show that they agree well with each other 
over the whole range of values of the scaled action S. 
The values of the local variance a^(w, Aw) can therefore 
be reproduced with an excellent precision by the leading 
order contribution alone. 

The expression of the spectral function C rr (w, Aw) in 
the semiclassical regime is obtained much in the same 
way as the expression of the spectral function D r (w), 
Eq. (p0|). Indeed, the contribution of a stable periodic 
orbit and all its repetitions to C rr (w) has an expression 
which is analogous to the one given in Eq. @. The only 
difference is that the amplitude A r is now replaced with 
the Fourier transform of the autocorrelation function of 
the Weyl transform of the perturbing operator along the 
stable periodic orbit |^|,^ 11 1. The formal resummation 
of all contributions leads to the semiclassical expression 
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FIG. 8. Same as Fig. 7 but in the case of the scaled energy 
E = -0.316. 

C rr (w,Aw) = ^MRe j°° dae^ Aw+i ^C Ar (a) , 

(30) 

with the autocorrelation function C A r (a) given by 



ds 



x (A{? + a/2) - (i) r ) (i(S' - a/2) - (i) r ) . (31) 

The autocorrelation function is a periodic function with 
the same period S r as the one of the orbit fllTTf. In or- 
der to check the validity of the semiclassical expression 
of C rr (w, Aw), it is easier to take the Fourier transform 
with respect to the scaled energy difference Aw on both 
sides of Eq. (|o|) and to integrate subsequently with re- 
spect to the scaled energy w. If TV is the number of reg- 
ular states in the used spectrum, one gets the following 
expression of the autocorrelation function 



C \r(S) = ^ E \(mr\A-(A) r \n r )f 



X COS 



(32) 



Fig. ^ compares the values of the autocorrelation func- 
tion, as computed with the help of Eq. (|32|) (full curve) 



and Eq. (f31j) (dashed curve). The comparison is done 
in the case E = —0.2 (N = 120) and over one period 
of the stable trajectory. It is seen that the agreement 
between both curves is good, except in the vicinity of 
the values S = and S — S r for which the expression 
in Eq. ([H]) is singular. The full curve would reproduce 
this singular behavior all the better as the number of 
regular states used in the computation would be larger. 
The finite number of regular states is also responsible for 
the observed little discrepancies between both curves in 
the considered interval of values of the scaled action. In 
spite of these discrepancies, the general good agreement 
between the full and the dashed curve allows to conclude 
that the semiclassical expression given in Eq. ( |30|) can 
be used to calculate the values of the spectral function 
C rr {w,Aw). By analogy with Eq. (|2^), one can also in- 
troduce a local variance cr^(w, Aw) corresponding to the 
distributions of non-diagonal matrix elements coupling 
regular eigenstates together, defined as 



a^(w, Aw) 



C rr (w, Aw) 

{poA w )) 2 



(33) 



It is to be noted that the spectral functions which are 
associated to the distributions of non-diagonal matrix el- 
ements involving chaotic (Eqs. ^?], |2^) and regular (Eqs. 
|30| , pl| ) states have similar classical expressions. As in 
the case of the diagonal matrix elements, each expres- 
sion takes into account the classical trajectory which is 
related to the component of phase-space connected to the 
studied subset of non-diagonal matrix elements. 

It has been checked previously that a very close esti- 
mate of the values of the total spectral function C (w, Aw) 
is obtained through the sum of the two spectral functions 
C cc (w, Aw) and C rr (w, Aw). Consequently, by virtue of 
Eqs. ^ and [53], the rescaled variance which is as- 
sociated to the distributions of all non-diagonal matrix 
elements can be written as 



a^(w,Aw)p 0it {w) 



/ — r cr c (w, Aw)po,c{ w ) 

poA w )J 

E^) a * {WjAw)poAw ). (34) 



As in the case of the mean value (A) t , Eq. (|22]), the 
contribution to the rescaled variance of (w, Awjpoj(w) 
of the ergodic (resp. regular) component of phase-space 
is weighted by the ratio of the mean density of states of 
the subset of chaotic (resp. regular) states to the mean 
total density of states. Equivalently, taking the Fourier 
transform on both sides of Eq. (|34| ) with respect to the 
scaled energy difference Aw, one obtains the following 
formula with the help of Eqs. |2TL |o| and |3| 



d(Aw) cos(AwS) cr t (w, Aw) po.t (w) 



P °' ciw) W(S) + (^)C A (S). ,sr„ 



poA w ) 



poA w ) 
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FIG. 9. Autocorrelation function C as a function of the 
scaled action S/2ir in the case of the scaled energy E — —0.2. 
Solid line: Computation with Eq. Dashed line: Computa- 
tion with Eq. 

The numerical comparison between both sides of this 
equation is done in Figs. 0b, c for the case E — —0.2. In 
Fig. |b, the full curve represents the Fourier transform 
of the rescaled variance whereas the dashed curve rep- 
resents the autocorrelation function C jj C (S) alone. One 
sees that both curves do not agree so much with each 
other. One has to compare Fig. with Fig. 0c in or- 
der to appreciate the quantitative improvement which is 
brought in the weighted sum by the contribution origi- 
nating from the stable periodic orbit. In this last figure, 
the full curve represents again the Fourier transform of 
the rescaled variance whereas the dashed curve repre- 
sents now the weighted sum of the autocorrelation func- 
tions on the r.h.s. of Eq. (fj5|). Contrary to Fig. 0b, both 
curves agree now well with each other over the whole 
range of values of S. The improvement is especially no- 
ticeable in the immediate vicinity of the peaked struc- 
tures located at positions which are multiple integers of 
S r /2n — 1.03305. This is due to the fact that the au- 
tocorrelation function C ' ^ r (S) contributes mostly to the 
Fourier transform in the immediate vicinity of these po- 
sitions, as shown by Fig. 0. Fig. 0c provides also an- 
other confirmation of the fact that the spectral functions 
C cr (w, Aw) and C rc (w, Aw) give really a negligible quan- 
titative contribution to the local variance. One can there- 
fore conclude that Eq. is able to reproduce the val- 
ues of the local variance with an excellent precision. This 
equation generalizes to the case of a scaled system with 
a mixed phase-space the expression of the local variance 
given in Refs. for a scaled system with an ergodic 
phase-space. As in the case of the diagonal matrix ele- 
ments, the generalization of Eq. (|34|) to a situation with 
several regular components in phase-space is straightfor- 
ward. Indeed, one has to add as many weighted rescaled 
variances a^(w, Aw)po tr (w) as there are stable periodic 
orbits in the underlying classical dynamics, as illustrated 
by the remaining figures. On the one hand, Fig. jsj^ (resp. 
Fig. [il]|a) compares the Fourier transform of the rescaled 
variance of (w, Aw)po t t(w) (full curve) with the autocor- 
relation function C ^ C (S) (dashed curve) in the case 




S/2n 

FIG. 10. Fourier transform of the rescaled variance 
of (w, Aw)po,t(w) of the distributions of non-diagonal tran- 
sition matrix elements as a function of the scaled action S/2n 
in the case of the scaled energy E = —0.4. (a) Solid line: 
Fourier transform of at (w, Aw)po,t(w); Dashed line: Auto- 
correlation function C ^ . (b) Solid line: Fourier transform of 
of (w, Aw)po,t(w); Dashed line: Weighted sum of autocorre- 
lation functions C ^ c and C ^ r . 

E = —0.316 (resp. E = —0.4). As in the previous case, 
there are many discrepancies between both curves be- 
cause the contributions coming from the stable periodic 
orbits are missing. On the other hand, Fig. |3|c (resp. 
Fig. |To|b) compares the Fourier transform of the rescaled 
variance (full curve) with the weighted sum of autocorre- 
lation functions (dashed curve) in the case E = —0.316 
(resp. E = —0.4). The sum now contains the properly 
weighted contributions originating from the stable peri- 
odic orbits which have been identified in the previous 
Section. In both cases, the very good agreement between 
the two curves confirms the validity of the generalization 
of Eq. (§4|). 

V. CONCLUSION 

This paper has been devoted to the study of the quan- 
titative contributions of the different components making 
up a mixed phase-space to the value of the mean charac- 
terizing the distribution of diagonal transition matrix el- 
ements and to the value of the variance characterizing the 
distributions of non-diagonal transition matrix elements. 
With the help of numerical computations done in the 
framework of the Hydrogen atom in a strong magnetic 
field, it has been shown that these contributions can be 
well identified in the semiclassical regime. The computa- 
tions have confirmed that the leading order contribution 
of the ergodic component to the mean is equal to the 
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average of the Weyl transform A of the perturbing oper- 
ator A along an arbitrary ergodic trajectory. They have 
also confirmed that the leading order contribution of the 
same component to the variance is proportional to the 
Fourier transform of the autocorrelation function of A, 
this autocorrelation function being also computed along 
an arbitrary ergodic trajectory. On the other hand, it has 
been found that the contribution of each regular compo- 
nent to the mean is equal to the average of A around the 
corresponding stable orbit. It has also been found that 
the contribution of each such component to the variance 
is proportional to the Fourier transform of the autocorre- 
lation function of A, this autocorrelation function being 
again computed around the corresponding stable orbit. 
For each studied quantity, the contributions coming from 
the ergodic and regular components have therefore sim- 
ilar expressions, each expression taking into account the 
particular classical trajectory which is related to the con- 
sidered component. The stable periodic orbits provide a 
convenient method to compute the contributions of the 
various regular components to mean and variance with 
high accuracy. This method is different from the one 
which has been proposed by Robnik and Prosen for the 
same purpose. As a final step, it has been shown that 
mean and variance can be expressed as a weighted sum 
of the contributions of all different components belong- 
ing to the mixed phase-space. The weight appearing in 
front of a given contribution has been identified as the 
ratio of the mean density of states of the corresponding 
component to the mean total density of states of the sys- 
tem. Although the study has been done for a particular 
scaling system, the results presented in this paper are rel- 
evant to all generic scaling systems with a small number 
of degrees of freedom having a mixed phase-space. 
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